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We consider the class of spin Hamiltonians on a ID chain with periodic boundary conditions 
that are (i) translational invariant, (ii) commuting and (iii) scale invariant, where by the latter we 
mean that the ground state degeneracy is independent of the system size. We correspond a directed 
graph to a Hamiltonian of this form and show that the structure of its ground space can be read 
from the cycles of the graph. We show that the ground state degeneracy is the only parameter 
that distinguishes the phases of these Hamiltonians. Our main tool in this paper is the idea of 
Bravyi and Vyalyi (2005) in using the representation theory of finite dimensional C*-aIgebras to 
^— ■») . study commuting Hamiltonians. 
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I. INTRODUCTION 



Classification of the phases of matter is a major problem at the heart of recent activities in quantum many- 

. body physics research, especially after the discovery of topologically ordered phases. The quantum double model of 

(— I ' Kitaev [1] and the string-net condensation model of Levin and Wen [2] show that the phase diagram of 2D systems can 
O |, be very complicated. These two models of commuting Hamiltonians are exactly solvable and exhibit a rich structure 
in their ground states as well as elementary excitations. In ID, however, we are not aware of any essential example 
! of ^ commuting Hamiltonian besides the Ising model. Moreover, simulation of the ground states based on tensor 
■ networks has been proven to be efficient in ID gapped systems [3]. So we expect to have an easier theory in ID. 
Before reviewing known results in this regard let us first describe what we exactly mean by a phase. 

Hamiltonians of interest are translational invariant spin Hamiltonians defined on a lattice with periodic boundary 
conditions. We assume that Hamiltonians are scale invariant, meaning that the ground state degeneracy is independent 
, of the system size, and also have a constant gap above the ground state energy which again is independent of the size 
' of the system. We say that two such Hamiltonians (their associated ground spaces) belong to the same phase if there 
exists a continuous path of Hamiltonians satisfying the above conditions starting with one of them and ending with 
the other. 

Chen et al. [4] and Schuch et al. [5] have recently shown that assuming the ground space of such a ID Hamiltonian is 
exactly describable by Matrix Product States (MPSs) [6] , the only parameter that distinguishes phases is the ground 
state degeneracy [12]. Although it seems that we can remove the assumption in this result by using the proven area 
law in ID [3], there are some technical obstacles. First, we are able to approximate the ground state of any gapped ID 
Hamiltonian by an MPS, whose bond dimension scales (polynomially) with the system size. In [4] and [5] however, the 
bond dimension is assumed to be fixed independent of the system size. Second, we need to show that the true ground 
state and the approximate MPS are in the same phase which again should be proved with a gap that is independent 
of the system size. 

?H ■ Yoshida [7] has considered the problem of classification of phases in the special case of stabilizer Hamiltoniairs. 
He has shown that any translational invariant and scale invariant stabilizer Hamiltonian in ID is equivalent to some 
independent copies of the Ising model, and in 2D we essentially have the 2D Ising model and the toric code. He 
has proved a similar result in 3D assuming some extra conditions [8]. Although stabilizer Hamiltonians seem very 
restrictive, they have beeir extensively studied in the theory of quantum error correcting codes aird quantum memories 
because of their simple structure of syirdrome measurements. The advantage of these results comparing to [4, 5] is 
that Yoshida's classification gives a characterization of logical operators, and consequently low anergy excitations as 
well. We emphasis on the latter because for instance in the models of quantum double and string-net condensation 
the structures of elementary excitations are much richer than that of ground states. 

In this paper we generalize Yoshida's results in ID. We consider frustration free commuting Hamiltonians that are 
both translational and scale invariant. Due to the commutativity assumption these Hamiltonians are automatically 
gapped. We associate a graph with a Hamiltonian of this form and show that cycles of this graph encodes ground 
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space of the Hamiltonian. We also comment that the low energy excited states are described by paths of the graph. 
Our main result is that the phases of the ground states of these systems are determined only by their degeneracy. 

Although our assumption of commutativity comparing to that of [4, 5] seems very restrictive, this at least can be 
easily verified. Even the scale invariance of a Hamiltonians is hard to check in general, but for commuting Hamiltonians 
we provide a simple criterion to verify that. Moreover, in the last section we show that our results can be applied to 
a wider class of Hamiltonians than commuting ones. 



Let us first fix some notations. Hilbert spaces are denoted by and Ha is the Hilbert space corresponding to 
register A. 1^ is the identity operator of T-La, and L(H^) is the space of linear operators acting on T-La- We may 
consider L(?^^) as a Hilbert space equipped with the Hilbert-Schmidt inner product. 

In this paper we are interested in the phases of ground spaces of Hamiltonians define on a ID chain with periodic 
boundary condisions. We say that ground spaces of two Hamiltonians belong to the same phase if there exists a 
continuous path of gapped Hamiltonians starting with one of them and ending with the other. Here we often consider 
translational and scale invariant Hamiltonians in which case we further assume that all Hamiltonians along the path 
are translational and scale invariant and their gap is independent of the system size. We do not impose this constraint 
for commuting Hamiltonians, i.e., two commuting Hamiltonians are said to belong to the same phase even if the path 
connecting them consists of non-commuting ones. This is because we are considering phases of ground spaces, and 
a phase is usually defined by translational and scale invariant Hamiltonians. So these two properties are essential in 
the definition of the problem but being commuting is an extra constraint which we impose for simplification. 

By this definition when two Hamiltonians H and H' have the same ground spaces, they belong to the same phase. 
In fact, the linear interpolation tH + (1 — t)H' , < t < 1, is the path between them and one can easily show that 
these Hamiltonians for all t are gapped. Local unitaries also do not change the phase of a ground space. That 
is Hamiltonians H = ^^^"^ ~ J2j ^^jj+i^'' ^ unitary U belong to the same phase, and the path 

H{t) = J2j c**'^^jj+ie~**'^, < t < 1, where U = e*'^, connects the two Hamiltonians. Note that translational and 
scale invariance are preserved in these two interpolations. We may also embed the Hilbert space of local spins into 
larger Hilbert spaces (and modify the local terms of the Hamiltonian). This does not change the ground space and 
we again obtain the same phase. 

Here we focus only on frustration free Hamiltonians. Thus by a shift of energy we may assume that the local 
terms of the Hamiltonian are positive semidefinite and the ground space energy is zero. Let H — be such 

a Hamiltonian where /ijj+i acts only on sites j and j + 1- Let Pjj+i be the orthogonal projection onto the positive 
eigenspaces of /ijj+i, i.e., Pj.j+i is a projection and kerhj j^i = kciPjj+i. It is straightforward to see that H is 
gapped if and only if H' = Pj.j+i is gapped. Moreover, H and H' have exactly the same ground spaces and by 
the above observation belong to the same phase. This means that without loss of generality we may assume that the 
local terms of the Hamiltonians are projections. Note that if H is commuting then H' is commuting as well. This is 
because Pj.j+i belongs to the algebra generated by hjj+i and the identity operator. 

The following lemma due to Bravyi and Vyalyi [9] is a key step towards understanding the ground space of commuting 
Hamiltonians. For the sake of completeness we also include a sketch of the proof. 

Lemma II. 1. Suppose that Xab and Ybc o,re two hermitian operators acting on finite dimensional Hilbert spaces 
T-La ®'Hb and T-Lb ® 'He, respectively. Assume that Xab ® Ac o,nd 1a ® Ybc commute. Then there exists an index 
set V and decomposition 



such that Xab o,cts nontrivially only on H^, and not Tip, (for all (3 Cz V), and similarly Ybc acts nontrivially only 
on Hfj,. More precisely, if we let Up be the orthogonal projection onto the subspace Hp^ ®'Hp, C Hb, then 



II. SOME KEY OBSERVATIONS 



Hb = ^np,®np,, 



(1) 



Xab 



Ybc 



i^A ® Up) Xab {1a ® Up) 

-.^iUp® Ic) Ybc {Up ® Ic) 

P 



^Rap, ® Ip,, 



Ip, ® Sp^c, 

P 



for some hermitian operators RaPi and Sp^c acting on Ha ® Hp^ and Hp, ® He respectively. 
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Sketch of proof. Using the Schmidt decomposition of Xab and Ybc as vectors in bipartite Hilbert spaces ^{Ha) "8) 
L('Hb) and L('Hb) ® L('Hc), and writing down the constraint [Xab ® Ic? 1a ^00] = we find that it suffices to 
prove the following statement. 

Let S and S' be two subsets of hermitian matrices acting on T-Lb such that all operators in S commute 
with elements of S' . Then there exists a decomposition (f) such that elements of S (S' ) act trivially on 

So we need to prove the above statement. Let Z{S) be the space of all matrices that commute with elements of S. 
Since elements of S are hermitian, Z{S) is a C*-algebra. and it contains S' . Using the representation theory of finite 
dimensional C*-algebras [10] we find that there exists a decomposition (1) such that 

Thus elements of S' have the desired form. Furthermore, it is not hard to see that a matrix commuting with such 
Z{S) has to act trivially onT-Lp^. □ 



III. GRAPHS ENCODE GROUND STATES OF ID COMMUTING HAMILTONIANS 

Consider a ID chain of N spins of finite dimension d (qudits) with periodic boundary conditions. We let the 
Hamiltonian be ff^v = Pj-j+i^ where Pjj+i acts only on sites j,j + 1 (Hereafter we adopt indices to be modulo 

N). By the observation of the previous section, since we let Hn to be frustration free, we further assume that Pjj+i's 
are projections. The extra assumptions that we impose in this section are that is translational invariant, and 
that Hn is commuting. 

According to Lemma ILl since Pj-i.j and Pj.j+i commute, there exists a decomposition 

C'-^H^,^H^. (2) 

of the Hilbert space of the j-th qudit such that Pj-ij {Pj.j+i) acts trivially on H^j {T~L^i). Using the fact that 
Pjj+iS are translations of each other, we may consider the same decomposition at site j + 1. As a result, Pj.j+i acts 
trivially on T-L^j+i as well. Therefore, Pjj+i can be written as 

where Q^j ^j+i acts on Ti^j T^^j+i and is a projection. Since Pj^+i's are shifts of each other we may ignore the 
index j and represent Q^j ^j+i by Qar,Pi- The combined Hilbert space of all TV qudits then can be written as 

TV 

Qi,...,Q«ev j=i 

We correspond a directed graph G to our Hamiltonian. The vertex set of G is the index set V of decomposition (2). 
A directed edge is drawn from a ^ V to /3 E V if ker Qar,Pi ^ 'H-a,. ®T~Lpi is non-zero. Note that this graph may contain 
loops, i.e., edges from a vertex to itself. Moreover, G is defined based on the local terms Pjj+i and is independent of 
the system size N . 

G encodes the ground space of the Hamiltonian as follows. Let (a^, . . . , a^) be a directed cycle of length in G 
which may contain a vertex or edge several times. Since — ?> a^'^^ is an edge of G, there is a non-zero vector \tfj) 
in kerQ^j ^j+i for j — 1, . . . , N. By decomposition (4), \ip) = \ipi) ® ■ ■ ■ ® I^n) is a state in the Hilbert space of N 
spins. Then using (3) we have 

HnIv) = ^Pjj+i\(pi) «) • • • (8) \(Pn) 
j 

= ^ Ifi) «) • • • (g) Q„3 \(Pj) «) • • • (8> \(Pn) 
j 

= 0. 
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Figure 1: An A''-qudit matrix product state (with periodic boundary conditions and) with bond dimension x is obtained by a 
chain of A*' maximaUy entangled states |<1>) of local dimension x followed by local maps Sj : ® C*. Then the final state 

equals {Si ® • • • <8> Sn) |$) (8> • • • ® |"1?). For translational invariant MPSs we have Si = S2 = ■ ■ ■ ~ Sn- From this representation 
it is clear that if \ip) is an MPS, then {Xi ® • • • (S)Xn) |<^) is also an MPS with the same bond dimension. Moreover, the 
maximally entangled states |$) can be replaced with any bipartite state with local dimension x smd we still obtain an MPS 
with bond dimension x- 



So \(pi) (g) • • • eg) \(Pn) is ill the ground space of H^. Indeed, all ground states of are of this form. From (3) and (4) 
we have 

N 

a\...,a"eV J = l 

kerQ^j ^i+i is non-zero iS a? o?^^ is an edge of G. Thus if we let Cat to be the set of ordered cycles of length TV 
in G, we have 

JV 

(Ql,...,Q«)eCjv j = i 

and 

AT 

dinikeriJAr= ^ JJ^ dimker Q^j . 

(Qi,...,a")eCN J = l 

Here we must consider ordered cycles because we have the freedom of choosing which vertex corresponds to the first 
site, so for example in C2, («,/?) and (/3,a) are counted as different cycles if a 7^ /3. 

As a summary we conclude that the ground space of -//at, for all A^, is completely characterized by the graph G and 
the subspaces kerQa^.^, corresponding to each edge a /3. Similarly to the ground states that correspond to cycles 
of G, it is easy to see that excited states with energy 1 correspond to paths of G, and eigenstates with eigenvalue 2 
correspond to unions of two paths. 

Note that the translational invariance assumption in this result is not crucial. If i?jv is not translational invariant, 
instead of a single index set V we have A^ index sets Vi, . . . , V^r; Vj corresponds to the j-th local term. Then the 
associanted graph would be A^-partite where there are directed edges from vertices in Vj to vertices in Vj+i. 

An implication of our result is that the ground states of commuting Hamiltonians can be represented by MPSs. An 
MPS representation consists of a chain of maximally entangled states on which we apply some local linear transfor- 
mations [11] (see Figure 1). To write the state jt/?) = \<~p\) \'{>n) defined above as an MPS, we need only to map 
the j-th maximally entangled state to \<^j). Thus \ip) is an MPS whose bond dimension is a constant independent of 
N. 

Let us work out an example to clarify the ideas. Let d = 4, and using = Cg) represent each site with two 
spin-half particles. Define 

where and az are Pauli matrices. Then — i^ commuting and translational invariant. The index set 

V associated with has four elements and components of decomposition (2) are 

(g7^„. =C(|00) + |11)), 

-«ft(g-H;3„=C(|01)-|10)), 

n-y.^n^^ = c(|oo) - 111)), 

H0^®He^ = <C{\Ql) + \lQ)). 



Figure 2: The graph associated with the local projection Pj,j+i = (Ij Ij+i ~ {<^x ® o-x)j ® (fz ® (jz)j+i)/2, and Hamiltonian 
ffjv — X]j ^jj+i- This Hamiltonian is commuting and translational invariant but not scale invariant. 



Note that in this case since these subspaces are all one dimensional, we do not need to specify subsystems "Han 
"Hq^,... individually. The corresponding graph G has four vertices and is depicted in Figure 2. For iV = 3 for 
instance, (a, 7,0) is a cycle of length three and then ([00) + |11)) ® (|00) — |11)) <E) (|01) + |10)) is a ground state. 
However, {a, 7, /?) is a path which cannot be completed to a cycle and then (|00) + |11)) (|00) — |11))®(|01) — |10)) 
is an excited state with energy 1. 



IV. SCALE INVARIANCE 



We now apply the assumption that the system is scale invariant, meaning that for sufficiently large N the ground 
state degeneracy is a constant independent of N. So we let Hjsf be commuting, and translational and scale invariant. 
Our goal in this section is to show that for every N 

dimkeriJAT — |Ci|, (5) 

where Ci is the set of loops of G. 

Let Di S Cn be a cycle of G. By going around Di many times, we may assume with no loss of generality that N, 
the length of this cycle, is sufficiently large. Let Cn = {Di, ■ • • , Dk} and Cn+i = {Ei, . . . , Ek'}. We have 

k 

dim ker Hn = JJ^ dimkerQo,^^^, , (6) 

i=l a^PeDi 
k' 

dimkeriJjv-i-i — JJ^ dimker Qa^./S, , (J) 

i=l a~^PeEi 

where by a — ?> /3 G we mean that a — J> /3 is an edge in the cycle Di, and similarly for a — > /3 G For a cycle D 
of length N let be the cycle of length tN which starting from the same vertex as D, goes t times around it. Then 
D^^^, . . . , Dj^~^^ are different cycles in Cn{n+i)j and similarly , . . . , Ej^, e Cjv(Ar+i) are distinct. Therefore, 

k 

dimkcriJAr(jv+i) > ^ Y[ {dimker Qa^^pi)'^'^^ , (8) 

1=1 a^PeDi 
k' 

dimkeri7Ar(jv+i) > ^ Y[ (dimkerQ^^^/j,)^ . (9) 

Now applying the scale invariance assumption we have dimker_ffjv(Ar+i) = dimkcriJAr = dimker_ffjv+i. So by 
comparing equations (6)-(9) we obtain 

dimkerQct^^^, = f (10) 

for every edge a /3 in one of the cycles Di or Ei. (Note that when a ^ /3 is an edge of G, dimkerQ^^^^, is 
non-zero.) We further conclude that k = k' ^ dimkcri?jv(Ar+i) and indeed 

Cn(n+i) = . . . , D^^^} = {E^, E^}. 
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Thus for some i, D^~^^ = . By the fact that N and A'' + 1 are relatively prime, we find that the arbitrarily chosen 
cycle Di consists of the repetition of a single loop. We conclude that the scale invariance assumption implies that all 
cycles of G are essentially loops. Then (5) is an immediate consequence of (10). 

The structure of the ground space of iJ^r in even simpler with the scale invariance assumption. Let q — > a be a 
loop of G, and let \'Pa) G T~La^ ^Hai be a vector that spans ker Q^^^q,, . Then \ipa)®^ is a ground state, and the 
ground space of Hm is spanned by these vectors for all loops a — > a of G. The vectors \tfa)'^^ are still MPSs and in 
fact translational invariant MPSs. 

The example of Figure 2 is not scale invariant because the graph contains cycles that do not come from loops. For 
the special case of the ID Ising model, the graph G consists of two vertices with two loops and no other edges. From 
these two loops and the above construction we obtain the all spin-up and all spin-down states as the ground states of 
the Ising model. 

Elementary excitations are also easier to characterize with the scale invariance assumption. As mentioned in the 
previous section, elementary excitations correspond to paths of G, and excitations with energy 2 come from unions 
of two paths. A sufficiently long path in a finite graph must contain a cycle. Since loops are the only cycles of such 
G, any elementary excitation of Hn, for sufficiently large N , contains several copies of states for loops a — J> a, 

in their subsystems. In fact an energy 1 eigenstate corresponds to a path of the form either a— >Q!---— J'/Jor 
/3 — > a a • • ■ — > a. For example, Ising model does not have any energy 1 eigenstate because its corresponding graph 
is a union of two loops and there is no edge of the form a — > /3 for a ^ p. 

V. PHASES ARE DISTINGUISHED BY THE GROUND STATE DEGENERACY 

We are now ready to study phases of ID spin chains with commuting Hamiltonians. Objects of interest are ground 
spaces of ID commuting Hamiltonians that are both translational and scale invariant. By the discussion of Section II 
without loss of generality we can also assume that the Hamiltonian is a summation of local projections. Thus we may 
use results of the previous two sections which are summarized as 

keri7Ar= (kerg,^,^,)^'^, (11) 

where kerQa^^a; is one dimensional and is spanned by \^a)- 

Here since we are interested only in the ground spaces we can replace the other Qar,i3i 's that do not appear in the 
above expression with the identity operator. More precisely, for every a ^ or a — (3 where a — >■ a is not a loop, 
define Qq,.,^, = la^ ® and for loops a — > a e Ci let Qa^^ai = Qq^.q,- Furthermore, define 

and Hn ~ '^j Pj,j+i- Then is still commuting and translational invariant. Moreover, its corresponding graph 
consists of the same loops as G, but no other edges. In fact, using (11), Hn has the same ground space as Hn and 
they are in the same phase. So we may replace Hn with Hn, i.e., we assume that Pjj+i — Pj,j+i and Hn = Hn- 

For every loop a — a fix arbitrary states |^ct,.) G Tia,^ and ) G T-tar Let Ua be a unitary operator that acts on 
TLa^ CS) Ha, in such a way that Ualifia) = ICctr) ® \^a,)- Recall that the Hilbert space of two qudits j and j + 1 can be 
decomposed as 

n^,+i n^^+i . (12) 

/3,76V' 

Thus we can define the two-qudit unitary operator U that acts as Ua on the subsystem/subspace H^i ^T-L^j+i when 
a — >■ a is a loop, and acts as the identity operator elsewhere. 

Define Pj j+i — UPj.j+iW , and H'j^ = Pj j^^. Since U is block-diagonal with respect to the decomposition (12) 
and acts trivially on and H j+i , the new Hamiltonian H'j^ is commuting. The corresponding graph of H'j^ is the 
same as that of Hn and consists of a union of loops. The only difference is that the state corresponding to the loop 
a — >■ a, is equal to U\(fa) = ICa^) <8) I^q;)- 

Two Hamiltonians Hn and H'j^ belong to the same phase because they differ only by local unitaries (see Section II) . 
So again without loss of generality we assume Pj.j+i — Pj j+i and Hn = H'n- In this step we turned I^^q), which 
belongs to Tia^ ® Ha, , into a product state. 
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By (11) the ground space of ff^v is spanned by vectors 

N N 

for loops a ^ a, where we set \a)j — |f j ) (g) IC^}- If we let TJ^v = J2i Pj-J+i such that 

then is commuting and translational invariant, and has the same ground space as Hff. Then they belong to 
the same phase. Moreover, since I^q^) and \^ai) were arbitrarily chosen, the only parameter that determines is 
the size of Ci, i.e.. the ground space degeneracy. For example for two loops we obtain the Ising model whose local 
projections are given by 

F = 1 ® 1 - (|0)(0| eg) |0)(0| + |1)(1| ® |1)(1|) . 

We conclude that the phases of ground spaces of translational and scale invariant commuting Hamiltonians in ID are 
characterized by their degeneracy. 

VI. SUMMARY AND OUTLOOK 

In this paper we described the structure of the ground states of translational and scale invariant, ID commuting 
Hamiltonians. We associated a graph with a commuting Hamiltonian which encodes the ground space in its cycles and 
the low energy states in its paths. Our results generalize Yoshida's work in ID who considers stabilizer Hamiltonians [7]. 
Comparing to [4, 5] instead of assuming that the ground states are described by MPSs, we imposed the assumption 
that the Hamiltonian is commuting. 

In Section IV we argued that ground states of the Hamiltonians under consideration can be exactly written as 
translational invariant MPSs. Thus we could have skipped the previous section and directly used the result of [4, 5] to 
conclude that the ground state degeneracy is the only parameter that distinguishes phases. Our arguments, however, 
are based on much simpler techniques and we preferred not to refer to [4, 5]. 

Our results can be applied on a larger class of Hamiltonians than the commuting ones. Let ~ be an 

arbitrary translational invariant (frustration free) Hamiltonian. As before without loss of generality we may assume 
that /ij j+i is positive semidefinite and the ground state energy of is zero. Suppose that there exists a positive 
definite matrix X such that 

{hj-ij (g) Ij+i) • (Ij-i g) Xj g) Ij+i) • (Ij-i gi hjj+i) = [tj-i (g hjj+i) ■ (Ij-i (g Xj g) Ij+i) • {hj-ij g) Ij+i). 

(13) 

Then the Hamiltonian H'^ = h'^ with local term 

/^L+i = i^]" ® X]l^h,,+,{X]I^ ® X]i^^ (14) 

is commuting. These local terms are still positive semidefinite, and one can easily observe that \ip) is a zero energy 
state of H'^ if and only if g) • ■ • g) x]^'^\lp) is the ground state of Hisi. Thus the new Hamiltonian is frustration 
free as well. Furthermore, since X is assumed to be positive definite and then invertible, the correspondence between 
ground states is one-to-one. Therefore, H'jq is scale invariant if and only if Hjq is scale invariant, and in this case 
results of our paper are applied. For instance, we showed that ground states of commuting Hamiltonians have MPS 
representations. Moreover, if is an MPS then X^^"^ ®- ■ ■(gjr]/^|(^) is an MPS as well (see Figure 1). As a summary, 
the existence of a positive definite matrix X satisfying (13) implies that the ground states of iJw = /ij j+i have 
MPS representations, and results of [4, 5] are applied. 

With the above technique sometimes we can turn a Hamiltonian whose ground space has MPS representation 
into a commuting one. In particular the converse of this observation holds for Hamiltonians with a unique injective 
MPS ground state in the following sense. Consider a translational invariant Hamiltonian with a unique translational 
invariant MPS ground state \(p). We assume that the bond dimension of this MPS representation is x ^-nd th^ 
corresponding map is 5* : g) C-^ — ?• C^. Then we have 

= (5' g) • • • g) S*) 1$) g) • • • g) 1$), 
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where |$) is the maximally entangled state with local dimension x- We may also replace S with US where f is a 
unitary map. In this case \ip) is replaced with U ® ■ ■ ■ ® U\ip) which by the discussion of Section II belongs to the 
same phase as \lp). 

S is injective because we assume that |(^) is an injective MPS. Then the map is well-defined on the support of 
S. In fact, we may identify the domain and support of S*, and by applying an appropriate U (replacing S with US) 
assume that S is hermitian and positive definite. 

We now introduce a commuting Hamiltonian as follows. Consider 2N spins of dimension x sitting on N sites of a 
chain. So there are two spins on each site which are denoted by I and r, and the Hilbert space corresponding to sites 
j and j + 1 is 

(Cf^.0C?,-)®(q,-+i)®C^Vi))- 

Define the local projection 

and let H = P^j+i. H is commuting, translational invariant, and frustration free with the unique ground state 
1$) (g) ■ • • ® |<I>). Moreover if we define 

- is-' ® s-^,)p,^,+ds-' ® s-^,i 

then H ~ has a unique ground state which is \ip). As a result, the Hamiltonian H with ground state \ip) 

can be turned into a commuting one using equation (14) with X^/^ = S. 

We conclude that MPS ground states and commuting Hamiltonians are in close relation via (14). The advantage 
of working with commuting Hamiltonians, however, is that verifying commutativity in general is much easier than 
checking whether the ground space is describable by MPSs. 

We leave the problem of classification of phases of ID systems for general Hamiltonians, without the commutativity 
assumption and that of [4, 5], for future works. 
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